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Abstract. Results are presented for some infinite series appearing in Feyn- 
man diagram calculations, many of which are similar to the Euler series. These 
include both one-, two- and three-dimensional series. The sums of these series 
can be evaluated with the help of various integral representations for hypcrgc- 
ometric functions, and expressed in terms of C(2), C(3); the Catalan constant 
G and Cl2(7r/3) where Cl2(S) is Clausen's function. 



1. Introduction 

High-precision calculations in Quantum Field Theory require the evaluation of 
Feynman amplitudes, which represent quantum corrections in a perturbative expan- 
sion. This expansion is governed by a small parameter, the fine-structure constant, 
and related to the number of loops in the associated Feynman diagrams. These are 
graphical devices to aid in the book-keeping. 

The evaluation of such Feynman amplitudes can technically be rather difficult, 
and one must often resort to approximate methods, valid in certain kinematical 
limits. For the case of electron-positron scattering, the full second-order corrections 
are not yet known. However, in the limit of high energies and small scattering 
angles, the problem simplifies: the relevant part of the so-called virtual amplitude 
can be expanded in powers of large logarithms. The coefficients and arguments of 
these logarithms can be extracted by the use of suitable Mellin transforms. We shall 
here review some of the techniques involved, and present some of the mathematical 
results of such calculations. 



2. Physics Problem 

The mathematical results we will discuss, derive from the "two-loop" Feynman 
diagram depicted in Fig. 1. The solid lines represent an electron and a positron 
scattering via the exchange of photons (wiggly lines). The momenta pi, p 2 , p[ and 
p' 2 only enter through the combinations 

(2.1) s=(p 1+ p 2 ) 2 , t={p 1 -p' l ) 2 = (p 2 -p' 2 ) 2 , m 2 =p 2 =p' 2 

where s and —t represent squares of the over-all energy and momentum transfer, 
respectively. Furthermore, m is the electron mass. The corresponding Feynman 
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pi — ^ — - — ^ — • — r— p'i 



P2 — — . — < — . — P ' 2 



Figure 1 . Double-box diagram 

amplitude is proportional to the integral 
(2.2) M= < dW(l-Ea*)A-»(a) 



+ crossed 



[D s {a)s + D t {a)t + D m (a)m 2 + D x (a)X 2 } 3 



where a fictitious photon mass, A, has been introduced in order to make the inte- 
grand well-defined for all aj. The D^a) and A(a) are homogeneous functions of 
seven a,-, one associated with each internal line in the diagram. 

Asymptotically, for s 3> \t\ ^> m 2 3> A 2 (high energies, small angles), the integral 
M has the structure 

(2.3) M ~ A 4 log 4 (- -0+^3 log 3 (- ••) + •■• 

where the arguments of the logarithms are given by the kinematical variables s, t, 
m 2 and A 2 . 

3. Mathematical Approach 

3.1. Mellin Transform. In order to extract, from the integral ([j~2]), the asymp- 
totic behaviour given by Eq. ( |2.3| ), we perform a Mellin transform. If the function 
f(x) behaves "almost" like x~ l at large x, then we define the Mellin transform as 

oo 

(3.1) J{z) = J x l ~ z f(x)dx 

o 

The fact that x l ~ z f(x) is only marginally integrable, is represented by poles in 
J(z), as z — > 0. The inverse transform is given by: 

c-Hoo 

(3.2) /(*) = _!__!_ / x *J(z)dz, c>0 

c— zoo 

where the contour of integration is as depicted in Fig. 2. 
It is instructive to consider the simple example 

(3.3) f(x)=6(x-l)J TI log k x, 
for which the Mellin transform gives 

( 3 - 4 ) AO = 

Thus, an expansion of the Mellin transform in inverse powers of £ is equivalent to 
an expansion in powers of logarithms. 
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Mellin contour 



-x — x — x — x- 



z plane 

Poles of integrand 

hi — x — x — x — x — x — x- 



Figure 2. Contour used for the Inverse Mellin transform 



3.2. Factorization Formula. In order to facilitate the integration over the cej 
in (2.2), we invoke the "factorization formula" [El |l2}| . For three terms, it can be 
written as: 

(A 1 + A 2 + A 3 )-p 

(35) [ — f dZ2 r ( Zl ) r ( Z2 ) r (P Z Zl Z z ^ a- Zi A~ Z2 A~ p+Zl+Z2 

1 ' ) J C1 2m J C2 2m T(p) 1 2 3 

Im Ai > 0, Cj > 0, p - ci - c 2 > 

where the integration contours labelled a,- run from Cj — zoo to cj + ioo. To prove 
the factorization formula, one will need two preliminary results. The first one is 

p — iiip/2 poo 

(3.6) u-P = -— / f- 1 e itu dt, Re p > 0, Im u > 0. 

r (p) io 

Proof: First, change variables by putting a; = —itu. Using the integral representa- 
tion of the Gamma function, we find, 

— / (-iuyPe^xP^dx = — e™P /2 u-PT(p) = U - p 

r (p) Jo r (w 

which is valid if u is purely imaginary with u/i > and Re p > 0. This result can 
be extended to any u with Im u > by analytic continuation. The second result 
needed is 

(3.7) e zT = — [ T(z)e l7rz/2 T- Z dz, c> 0, Im T > 0. 
Proof: Consider the Mellin transform of e lTx [ci. Eq. (3.6)], 

poo 

I' = x z - 1 e lTx dx = T{z)e ilTZ l 2 T- z , 

Jo 

which is valid for Re z > and Im T > 0. By taking the inverse Mellin transform 
of this expression, one should get e lTx , 



e iTx _ 

2m 



— ( x- z T{z)e l7 ' z ' 2 T- z dz. 
m J c , 
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By putting x = 1 in this expression, the desired result is obtained. One is now 
capable of proving the factorization formula. 
Proof of the factorization formula: 

(Ax +A 2 + A 3 )-p = —— / dttP~ 1 e^ Al+A2+A3 '> 

r 0) Jo 

—iirp/2 rOO 

&tt p - l e itM e itM e UA3 



r(p) 

r(p) io 7 C1 2™ i C2 2m 1 1; 1 2j 

x e" zi /2 e™ 22 /2 (t Ai ) ~ Z1 (tA 2 ) ~ 22 e u M 

xA- Zl A- Z2 T{p- z\ - z 2 y^ p ~ z i- z ^/ 2 A- p+z ^ 



, C1 ~„ Jcj 

xe l7rzl/2 e" 22/2 A7 z Mr Z2 / drt^ 1 ^ 1 " 22 ^ 3 



from which the factorization formula follows. The formula is valid for Im A t > 0, 
Cj > and p — ci — c 2 > 0. It is trivially generalized to (Ai + A 2 + • • • A„)~ p . 

3.3. Emergence of Sums. Performing a Mellin transform in s, with |i|=l, m 2 = 
s' ,m , A 2 = s nx , and r/\ < rj m < 0, one is led to integrals over z, of the kind 
given by Eq. fl3.5p . These integrals are most conveniently performed by the use of 
Cauchy's theorem, closing the contour in the left or right half-plane. Expanding 
these summands in powers of the Mellin transform variable, one is led to sums 
involving, as residues, the T function and its derivatives, 

4>(z) = (d/dz) io g r(z) = r'(2)/r(z) 

and 

i>'{z) = (d/dz)ip(z). 

The combination 7 + tp(n) frequently appears in the summand, and several series 
of this type are evaluated using the following integral representation, 

f 1 \-t z ~ Y 
7 + 'ijj(z) = dt- 



1-t 

Here, 7 = 0.577 216 ... is Euler's constant. 

3.4. Special Functions. Our evaluation of the sums makes use of several special 
functions. We shall list a few of them along with some of their properties here: 
Riemann's zeta function: 

The sums of several series can be expressed in terms of the Riemann £ function: 

00 

k=l 
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with 

G 



00 i 2 
^ p = ((2)4 = 1.644 934. 



OQ 

Ep = C(3) = 1-202 057.... 

fc=i 

Polylogarithms: 

Furthermore, we need Euler's polylogarithms: 

00 £ 

Li„(z) = ^^, \z\<l, n = 0,1,2,..., 
These functions also have integral representations: 

(, 8) u M = tS^[\ t >°Z"- 2 ^ 1 - t °\ „ = 2,3,4,... 
(n-2)! 7 i 

(3.9) =^ df Ll "- l(t) , n=l,2,3,..., 

with Li„(l) - {(n), n = 2, 3, 4, . . . , Li n (-1) - (1/2"- 1 - l)C(n), n = 2, 3, 4, . . . . 
Clausen's function: 

On the unit circle, the imaginary part of the dilogarithm is Clausen's function: 

Cl 2 W=Im[Li 2 (e i9 )]=f;^. 

fe=i 

Clausen's function has its maximum value for 8 = ir/3, Cl2(7r/3) = 1.014 942 
Furthermore, Cl2(7r/2) = G = 0.915 966 . . . , where G is Catalan's constant. 

3.5. A Theorem For Summing Series. Many of the results given here were 
obtained using "Theorem 1" of Q : 

(3-io) E^[7 + ^(i + M] = (y + |)c(3) + /l>i 2 (^) 

71 = 1 ^ ' j=l ^ ' 

n=l v 7 j=l v 

for k = 1, 2, 3, . . . . The sums over j are understood to be zero when k = 1. 

To prove the theorem, we use properties of Nielsen's generalized polylogarithms 
|6| . For a detailed proof of the theorem, we refer to |l(| . 



4. One-dimensional series 



We present a collection of one-dimensional series and their sums. For some of 
the series we also present the methods used to sum them. For the rest of the series, 
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the proofs can be found in 0, 0. 



oo 1 

(4.1) Yl T [7 + ^(1 +n)]=2C(3) 

n—l 

oo 

(4.2) ^^[ 7 + ^(l + 2n)] = T C(3) 

77=1 

oo „ 

(4.3) E + 4>(1 + 3n)] = 5C(3) - -£ci 2 (| 

n—l 

(4.4) J2 -[7 + ^(l + 4n)]--C(3)-27rG 

n—l 

(4-5) £^[7 + ^(1+ 6n)] = ^C(3) - ^Cl 2 (J 



n=l 



Some of these series were known already by Euler! Alternating series give similar 
results: 



( 46 ) E L ^ L [7 + V'(l+n)] = - 5 C(3) 



71 = 1 



n 2 L ' rv /J 8 



(I'm oo 

(4.7) £ + ^(1 + 2n)] = -C(3) - ttG 

n A 16 

71=1 

(4.8) E K -Jrh + Ml + 3n)] = -C(3) - 2^C1 2 (- 



77=1 



The results presented so far, all follow from the theorem or the immediate corollary 
thereof. 

Some series encountered had to be summed by other means: 



(4-9) E^ny[7+^a+^=c(2) 

77=1 ^ ' 

OO 

(4-10) E ra 2 (n + 1) [7 + ^(1 + »)] ^ 2C(3) - C(2) 

n—l ^ ' 
oo ^ 

(4.H) E /■ , i ^ [7 + ^(1+ n)]=-C(3) + C(2) 

77=1 X 7 



Related sums, where the argument of the psi function is shifted by a small integer, 

are easily obtained using the recursion relation ^(1 + n) = ^->{n) + 1/n. 

Proof of $4.!\): We apply the integral representation of the psi function before 
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summing over n: 



7i—l ■ ' n—1 x ' u 



dt- 



l - 1 



dt- 



1 



1 - 1 



dt- 



1 



1 - 1 



1- 1 2^(1,1; 3; t) 



l-i[t + (l-t) log(l-t)] 



= C(2) 



Here, we used (7.3.2.150) of Prudnikov pj to rewrite 2-Fi- The resulting integral 
is just the definition of Li2(l) which is well known. 

Series containing the derivative of the psi function were also encountered: 



(4.12) 
(4.13) 
(4.14) 
(4.15) 



£-V'(n)=2C(3) 

n=l 

oo 

£-V'(l+n)=C(3) 



oo 1 

^-fn(n+l) v v ' 

oo 1 
n—1 v y 



Proof of Series (4-.lt ): We start by using the integral representation of the trigamma 
function: 



logf 



oo 1 oo „i 

n=l n=l 1/0 

= ^ ^logilog(l-t) 



■logtlog(l-t) = C(3) 



The resulting integral is just the integral representation of Lis(l) which equals C(3). 
Also some series bilinear in ip(z) can be summed by similar methods: 

OO 

^n(n + l) 

oo 1 

( 4 ' 1? ) E 7 T~TT b + ^)}h + ^(1 + n)] = C(3) + C(2) 

^ «(n + l) 

oo 

( 4 - 18 ) E n(n + n b + VX")][7 + ^(2 + n)] = 3 
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(4.19) 
(4.20) 
(4.21) 



oo 1 

E^TI) b + M + »)] a = 3C(3) 

71—1 V / 



-[ 7 + i>{! + n)][ 7 + ^(2 + »)] = 2C(3) + C(2) 



V : 

^ n(n + 1) 

71=1 v ' 
oo 1 

E^n)[7 + ^(2 + ,)f = c(2) + 3 

n— 1 v / 



These series are summed by using the integral representation of each of the factors 
7 + ip(l + n) before summing over n. The resulting two-dimensional integral yields 
the desired result. 

Some one-dimensional series containing T-functions also appear. Consider the 
following series: 



(4.22) 



(4.23) 



ylM = __v (3) 

^ n 2 r(2n) 3^ ; 

n— 1 v ' 



4?r 



<*(§) 



T(2n) 



Proof of (4.2i): 



V 1 F(n| _, v 
^ 7i 2 r(2n) ^n 2 r(i+n) 

n— 1 v y n— 1 / 



i r(|)r(n) /i 



,r. j I. i. 1.1; ?,2,2; \ 



dl .F, j 1,1,1; I 2; it 



Next, we make use of (7.4.2.353) of Il4l and the substitution yftjl = sin ¥ to get 



x di 



• 2 >/* o 

— arcsm — = 2 
i 2 



du 



2 tan | 



Integration by parts yields: 



2 m 2 log ( 2 sin - 



du u lot 



du it log 



(2 sin I) 



( 2sin ?) 



,CM + ^C1 2 (| 



3 ,v ' 3 

after using (6.46) and correcting the misprint in (6.52) of ||. 

5. Two-dimensional series 

In the evaluation of Feynman diagrams by the above outlined methods, also two- 
and three-dimensional series are required. For the purpose of summing certain two- 
dimensional series, the following results are very useful: 



(5.1) 



oo 1 

E — 

^ n 4- 



1 



1 



n + a n + b b — a 



+ b) - -0(1 + a)] , a^b 
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and 

(52) g T(n + k) _ T(k) 



n=l 



T(l + n + 2k) r(l+2fe) 



Some selected two-dimensional series that do not involve the psi function are: 

oo oo 1 

( 5 - 3 ) yy^i — rr = 2 c(3) 

^ ^ nk(n + k) w 



n=l fc=l 



x x . 

<"> g o g M „ +fe)(1+ „ +fc) - 



OO OO 



71=1 fe=l V ' ' 

(5 6) i r(2fc)r(n + fc) _ 

ra=lfc=l v ; 
A few typical series involving the psi function are: 

(5J) ^ 2 /fc(nVfcKl + n + jfe) = 2C(3) 

n=0 fc=l v / v / 

(5 101 V V 7 + ^(l + n + 2fc) _ 7 

6. Summary 

We have seen that Quantum Field Theory relates to rather interesting mathe- 
matics. In physical problems of central interest (see, e.g., |I| ||, [| [l3|, |l5)), series 
related to the Euler series emerge. They have been summed in terms of C(2), C(3) 
and closely related irrational constants. A computerized treatment of such sums is 
also available |l6| ]. 

It is interesting to note that the transcendental basis that is required for the 
considered types of sums is rather restricted. In particular, there are no £(4), C(5), 
etc. The restrictions of such transcendental bases have received a lot of attention 
recently, in particular by Broadhurst and Kreimer (|, fjj (see also ||), who relate 
these restrictions to properties of Hopf algebras. 
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